AD  Mo. 

DDC  FILE  COPY  AD  A 0 5 49  3 4 


V' 


FOR  FURTHfR  TRAH  ^ 


h ' NRL  Memorandum  Report  3641 


Effects  of  Collisions  on 
Level  Populations  and 
Dielectronic  Recombination  Rates 
of  Multiply-Charged  Ions 

V.  L.  Jacobs  and  J.  Davis 
Plasma  Physics  Division 


NAVAL  RESEARCH  LABORATORY 
Washington,  D.C. 


Appruved  for  public  release:  distribution  unlimited. 


sECuftiTv  Classification  of  this  page  ^Pii*n  Oai>  Enw^) 


REPORT  DOCUMENTATION  PAGE 


T BEPOST  number 

NRL  Memorandum  Report  3641 


12.  OOVT  ACCESSION  NO. 


^ ^ title  f«n<l  : 

f u > IjFfects 
y^D  ME 
t^iJLWL 


title  m— — 

FECTS  OF  ELUSIONS  ON 
CTR( 


1C  MC 


^VEL  mPULATIONS 
IlNATIO^ 


J3ATES  OF 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3 RECIPIENT'S  CATALOS  NUMBER 


Interim  repMti 


IPERIOD  COVERED 

: a continuing 


E.  PERPORMING  ORS.  REPORT  NUMBER 
»■  CONTRACT  OR  ORAy  WUMimH, 


9.  PERTORMINS  ORGANIZATION  NAME  AND  ADDRESS 


Naval  Research  Laboratory 
Washington,  D.  C.  20375 


1.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

National  Aeronautics  and  Space  Administration 
Washington,  D.  C.  20545 

T5  MONITCPING  AGENCY  NAME  A ’AOOPESS^//  difUr^t  fi 


TT-wraairTg  TCfSSHTTptSjW^TTBK 

APKA  * UNIT  NUMtCNS 

NRL  Problem  No.  77H02-24 
NASA  Grant  No.  DPR-60404-G 


GES 

is:  SECUPITV  class,  for  ihi*  rmport: 

UNCLASSIFIED 


OeCLASSi  FI  CATION'  downgrading 
schedule 


16  DiSTRIOoTiON  STATEMENT  (ci  thit  Rmpon) 


Approved  for  public  release;  distribution  unlimited. 


KTWiailTiftM 


MmA  20.  U fro«i  Ri 


IS-  SUPPLEMENTARY  NOTES 


19  KEY  WORDS  ^Conrlnu*  on  povoroo  oido  if  nocooooo'  identity  by  bfecA  niatboO 

Dielectronic  Recombination 
Ionization  Equilibrium 
Atomic  Processes  in  Plasmas 
Autoionization 


^20.  nABSTRACT  (Conttw  on  rovorao  mid*  if  nocoooAry  aid  idonNfy  by  blocb  maibor) 

A generalization  of  previously  reported  statistical  theories  is  developed  for  determining  the 
excited  level  populations  and  the  ionization-recombination  balance  of  multiply-charged  atomic 
ions  in  an  optically-thin  high-temperature  plasma.  Account  is  taken  of  the  most  important 
collisional  and  radiative  processes  involving  bound  anc\^|utojpjizing  levels  in  three  consecutive 
ionization  stages,  ^e-nlitoia  a set  of  rate  equations ^for  the  population  densities  of  the  low- 
lying  levels  which  contains  effective  excitation,  ionization,  and  recombination  rates  describing 


nn 

VU  I JAN  7J 


1473 


EDITION  OP  ■ NOV  AS  IS  OBSOLETE 
S^N  0I02-LF-0)4-6601 


SECURITY  CLASSIPICATIOH  OP  TniS  PAGE  fBBMi  DM.  ^IWPd) 


jA 


4 


SECURITY  CLASSiriCATION  OF  THIS  PAGE  rWAwi  0««  enl«r«<0 


20.  Abstract  (Cont’d) 

indirect  transitions  through  the  more  highly-excited  bound  and  autoionizing  levels.  The  familiar 
corona  model  equations  for  the  ground  state  populations  are  recovered  by  making  the  assumption 
that  all  excited  states  decay  by  only  spontaneous  radiative  or  autoionization  processes.  When 
collisional  processes  become  efficient  in  depopulating  the  highly-excited  levels  important  in 
dielectronic  recombination,  the  effective  rate  of  recombination  must  be  described  by  a coUisional- 
dielectronic  recombination  coefficient.  Results  of  calculations  are  presented  for  the  coUisional- 
dielectronic  recombination  rate  coefficients  for  recombination  of  F^**^  ions.  At, an 

electron  density  of  10^^  cm~^,  dielectromc  recombination  is  still  theaominant  recombination 
process.  However,  the  coUisional-dielectr^c  recombination  rate  coefficients  are  found  to  be 
reduced  by  about  an  order  of  magnitude  fiom  their  corona  model  values  due  to  the  effects  of 
multiple  collisional  excitations  on  the  populations  of  the  highly-excited  bound  levels  of  the 
combined  ion.  The  dielectronic  recombination  rates  into  these  highly-excited  levels  are  found 
to  be  enhanced  by  the  effects  of  collisionally induced  angular  momentum  redistribution  on 
the  populations  of  the  autoionizing  levels.  \ 
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EFFECTS  OF  COLLISIONS  ON  LEVELS  POPULATIONS 
AND  DIELECTRONIC  RECOMBINATION  RATES 
OF  MULTIPLY  CHARGED  IONS 


I.  INTRODUCTION 

The  physical  properties  of  high-temperature  laboratory  and  astrophysical  plasmas  can  be 
deduced  from  an  analysis  of  optically-thin  spectral  line  intensities  emitted  by  multiply-charged 
atomic  ions  undergoing  transitions  between  low-lying  bound  levels.*’^  The  theoretical  predic- 
tion of  the  spectral  line  intensities  is  complicated  by  the  fact  that  the  low-lying  levels  can  be 
significantly  populated  by  cascade  transitions  from  the  more  highly-excited  levels.  In  the  solar 
corona  and  in  low-density  laboratory  plasmas  such  as  tokamak  discharges,  incompletely- 
stripped  ions  can  recombine  predominantly  by  means  of  the  dielectronic  recombination  pro- 
cess described  by  Burgess^,  in  which  highly-excited  bound  and  autoionizing  levels  can  play  the 
most  important  role.  In  this  investigation,  a consistent  statistical  description  is  given  of  the 
contributions  from  highly-excited  bound  and  autoionizing  levels  to  effective  excitation,  ioniza- 
tion, and  recombination  rates  between  low-lying  bound  levels.  Calculations  are  presented  for 
recombination  rates  for  which  these  contributions  are  expected  to  be  most  important. 

The  dielectronic  recombination  process  described  by  Burgess  ^ is  initiated  with  the  forma- 
tion of  a doubly-excited  autoionizing  state  through  a radiationless  capture 


jK'  + tz+l)  (,)  + 


of  an  electron  into  a state  nl  accompanied  by  the  excitation  / — * / of  one  of  the  bound  elec- 
trons in  the  recombining  ion.  The  initial  state  i of  the  (z  -I-  l)-times  ionized  ion  js 


usually  the  ground  state  but  recombination  from  a metastable  state  may  also  be  important. 
Recombination  is  accomplished  when  the  doubly-excited  state  undergoes  a stabilizing  radiative 
transition 


X + (/'.  nD-^X-^  Ho.  fil)  +fho 

to  a singly-excited  state  below  the  ionization  threshold. 


(2) 


At  sufficiently  low  densities  all  stabilizing  radiative  transitions  (2)  will  be  followed  by  cas- 
cade transitions 


+ (/„,  n/)  - + (/„,«"/'')  +1>0," 


(3) 

which  will  eventually  terminate  in  the  ground  state  jg  — tiglg  or  a metastable  state  of  the 
recombined  ion.  The  total  recombination  rate  to  the  ground  state  is  then  given  simply  by 


Rd  - L >0-  «/> 

f ,nl 


(4) 


where  A,(/',  nl  — ig,  nl)  denotes  the  spontaneous  radiative  decay  rate  for  the  stabilizing  transi- 
tion (2).  The  doubly-excited  level  population  densities  nl)  are  given  by 

MO)  KC,an(' 

- n n'  ’ 

Agit , nl)  -t-  Ar(i , nl)  ,c\ 


where  C„p(/  — * ; , nl)  is  the  rate  coefficient  for  the  radiationless  capture  (1)  and  the  total  decay 
rates  due  to  all  spontaneous  radiative  and  autoionization  processes  are  denoted  by  A^d',  nl) 
and  AgO',  nl),  respectively.  and  M(i)  are  the  electron  and  the  initial  ion  densities,  respec- 
tively. 


It  has  been  pointed  out  by  Burgess  and  Summers^  that  for  sufficiently  high  electron  den- 
sities the  most  important  singly-excited  levels  will  be  significantly  depopulated  by  the  collision- 
excitation  and  ionization  processes 
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+ (v  «/)  + 6 - - + (/„,  nT)  + e- 

and 


(6) 
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(7) 


Bates,  Kingston,  and  McWhirter^  have  given  a definition  of  the  effective  recombination  rate  to 
the  ground  state,  for  which  they  proposed  the  name  "collisional-radiative  recombination  rate", 
when  the  singly-excited  level  population  densities  Niig,  nl)  are  determined  by  both  collisional 
and  radiative  processes.  However,  dielectronic  recombination  was  not  taken  into  account  in 
their  investigation.  Burgess  and  Summers^  proposed  the  name  "collisional-dielectronic  recom- 
bination rate"  for  the  effective  recombination  rate  obtained  with  the  inclusion  of  dielectronic 
recombination.  The  collisional-dielectronic  recombination  rate  is  given,  in  terms  of  the  total 
collisional-radiative  transition  rates 


“ L ^^>0’  'o-  "o'o'f  ^('o-  "/)• 

»'  (8) 

The  collisional  processes  (6)  and  (7)  together  with  their  inverses  tend  to  bring  the  singly- 

excited  level  population  densities  NO^,  nf)  closer  to  their  Boltzmann-Saha  equilibrium  values, 

which  are  found  to  be  much  smaller  than  predicted  when  only  radiative  and  dielectronic 

recombination  and  cascade  transitions  are  taken  into  account.  Therefore,  collisional  processes 

tend  to  reduce  the  effective  recombination  rate.  Since  the  most  important  An  ^ 0 collisional 

processes  are  found  to  be  associated  with  the  n'  n ± 1 transitions,  the  reduction  may  be  said 

to  be  accomplished  through  multiple  collisional  excitations  to  more  highly-excited  levels  from 

which  ionization  is  more  probable. 


Electron-  and  proton-induced  transitions  corresponding  to  — n and  /'  — / ± 1 are 
known  to  be  the  most  rapid  collisional  processes  involving  the  multiply  charged  ions  in  a 
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predominantly  hydrogen  plasma.  Burgess  and  Summers^  made  the  assumption  that  these 
processes  are  rapid  enough  to  establish  a Boltzmann  distribution  with  respect  to  / for  all  values 
of  n.  They  also  pointed  out  that  the  corresponding  angular  momentum  redistribution  processes 
involving  doubly-excited  states 

+ (/',  nl)  + e + (/',  n.l±l)+e- 

(9) 

X + (/',  nl)  -I-  //  ■*•  — A'  + O'.  n.t±  I ) + //  + 

(10) 

would  produce  an  enhancement  of  the  dielectronic  recombination  rate  into  the  sublevels  associ- 
ated with  a given  value  of  n,  but  this  effect  was  neglected  in  their  calculations.  Independent 
calculations  by  Jacobs,  Davis,  and  Kepple^  and  by  Seaton  and  Storey^  have  now  established  that 
collisionally-induced  angular  momentum  redistribution  of  the  doubly-excited  level  populations 
can  amplify  the  dielectronic  recombination  rates  into  the  highly-excited  levels.  The  calculations 
of  the  colisional-dielectronic  recombination  rates  reported  in  this  paper  are,  to  the  best  of  our 
knowledge,  the  first  in  which  explicit  account  is  taken  of  the  rates  for  collisionally-induced  an- 
gular momentum  mixing  of  both  the  singly-  and  the  doubly-excited  states. 

A detailed  discussion  is  given  in  section  II  of  the  individual  collisional  and  radiative 
processes  which  are  taken  into  account  in  the  present  investigation.  In  section  III  a set  of  rate 
equations  for  the  population  densities  of  the  low-lying  levels  is  derived  which  contains  effective 
transition  rates  describing  indirect  transitions  through  the  more  highly-excited  bound  and  au- 
toionizing  levels.  By  taking  into  account  processes  involving  bound  and  autoionizing  states  in 
three  consecutive  ionization  stages,  we  obtain  a generalization  of  previously  reported  statistical 
theories. Certain  aspects  of  the  general  theory  are  illustrated  in  section  IV,  where  calcula- 
tions are  described  for  the  collisional-dielectronic  recombination  rate  coefficients  for  recombi- 
nation of  Fe  — Fe  ions.  Our  conclusions  are  given  in  section  V. 


II.  INDIVIDUAL  COLLISIONAL  AND  RADIATIVE  PROCESSES 


In  this  section  we  discuss  the  individual  collisional  and  radiative  processes  which  are  ex- 
pected to  be  most  important  in  determining  the  level  populations  of  multiply-charged  ions.  In 
some  cases  we  introduce  total  transition  rates  or  rate  coefficients  between  levels  which  include 
the  contributions  from  both  collisional  and  radiative  processes.  Since  we  will  ignore  collisions 
between  the  multiply-charged  ions,  we  shall  assume  that  they  are  present  in  sufficiently  small 
concentrations  in  a predominantly  hydrogen  plasma.  Since  only  spontaneous  radiative 
processes  will  be  taken  into  account  in  the  present  investigation,  we  shall  also  assume  that  the 
plasma  is  optically-thin  to  its  own  emitted  radiation  and  that  no  external  radiation  fields  are 
present.  However,  we  will  indicate  which  of  the  total  transition  rates  and  rate  coefficients  can 
be  influenced  by  contributions  from  processes  induced  by  a radiation  held. 

In  the  present  investigation,  we  will  take  into  account  only  transitions  between  levels 
within  an  ionization  stage  and  those  connecting  adjacent  ionization  stages.  Double  ionization 
processes,  for  example,  will  be  neglected.  To  simplify  the  notation  used  in  the  following  sec- 
tion, we  will  adopt  the  convention  of  denoting  bound  states  of  the  three  consecutive  ions 
jC  + (z  + i)^  ^ + (z)^  gjjjj  j + (z-l)  jjy  indices  /,  7,  and  k,  respectively.  Autoionizing  states  of 
the  ions  X and  ^ -1 ) jjg  denoted  by  a and  6,  respectively.  Different  states  belong- 

ing to  the  same  stage  of  ionization  will  be  distinguished  by  using,  for  example,  7, 7',  and  7.  In 
this  investigation  the  singly-excited  bound  states  will  be  specified  by  7 = /„,  n/,  where  ig  is  the 
ground  state  of  the  ion  \ The  corresponding  doubly-excited  autoionizing  states  will  be 

specified  by  a — /',  nl,  where  /'  is  an  excited  state. 

The  total  collisional-radiative  transition  rate  connecting  the  bound  states  7 and  7',  which 
are  in  the  same  ionization  stage,  is  given  by 

A.iJ'^J)  0(Ej.  -Ej)  + C^(7'-7)  + C^Aj'-j). 
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The  rate  for  the  spontaneous  radiative  decay  process 

x + + (y)  +ii<o 


(12) 


is  denoted  by  A^U'--  y),  and  the  9 function  is  used  to  insure  that  emission  terms  contribute 
only  when  Ej  < Ej.  The  corresponding  stimulated  emission  and  photo-excitation  rates  must 
be  included  when  the  interaction  with  a radiation  field  is  investigated.  The  rate  coefficients 
describing  the  collisional  processes 

X + (^Hj')  + 6 --  + (y)  -I-  e ■ 

(13) 


and 


T + (y')  -I- //+- + (y)  -I-//  + 


(14) 

which  are  denoted  by  CgiJ'-^  J)  and  C^+  (y'  — y),  are  functions  of  the  local  temperatures 
and  T^  + . The  electron  and  proton  densities  are  denoted  by  and  + . Electrons  are 
known  to  be  more  efficient  in  causing  transitions  corresponding  to  energy  differences  compar- 
able with  kgTg,  while  protons  become  more  important  for  transitions  between  the  nearly- 
degenerate  /-sublevels  corresponding  to  large  w-values.  With  increasing  n,  radiative  decay  be- 
comes less  probable,  whereas  collisional  transitions  between  neighboring  levels  occur  more  ra- 
pidly. 


The  bound-free  and  total  free-bound  rate  coefficients 


5(/,y)  “S^fy-/) 


(15) 


and 

a (y,  ;)  •=  a (/  — y)  -I-  03  (/  — * y ) 

(16) 

describe  direct  transitions  between  the  bound  states  ; and  y,  which  are  in  adjacent  stages  of 
ionization.  The  rate  coefficients  for  the  electron  impact  ionization  and  three-body  recombina- 
tion processes 
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jf  + u)(y)  + f -z:  + 


(/)  +e  + e 


(17) 

are  denoted  by  S,  (y  — » /)  and  03  (/  -*y),  respectively.  Charge  exchange  processes  are  neglect- 
ed in  the  present  investigation.  The  rate  coefficient  describing  the  spontaneous  radiative 
recombination  process 

A"*- <-+'>(/)  + e " — + (y  ) + if  <0 

(18) 

IS  denoted  by  a,(/— y).  To  investigate  the  interaction  with  a radiation  field,  S(i,j)  must  in- 
clude a term  describing  photoionization,  while  the  stimulated  radiative  recombination  rate 
coefficient  must  be  included  in  a (y,  /).  The  direct  ionization  and  radiative  recombination 
proce.'ses  result  predominantly  in  the  production  of  relatively  low-lying  states.  Three-body 
recombination  becomes  more  important  than  direct  radiative  recombination  with  increasing  n 
and  Ng,  but  dielectronic  recombination  is  more  often  the  dominant  recombination  process  at 
low  and  intermediate  densities. 


We  will  assume  that  transitions  between  the  autoionizing  state  a and  the  adjacent  contin- 
uum state  /,  e,  occur  only  by  means  of  the  radiationless  capture  and  inverse  autoionization 
processes 


j.  + (r-n)  (,)  + ^ -(€,)  = + (fl). 

The  capture  coefficient  and  the  autoionization  rate  will  be  denoted,  respectively,  by 


(19) 


and 


Cear^o-  ')  = C,a^(/  —a) 


(20) 


(/,  a)  = (a  -*  /). 

(21) 

We  will  neglect  the  corresponding  three-body  capture  and  collisional  ionization  processes 
described  by  Weisheit^,  which  may  become  important  only  at  very  high  densities  or  for  very 
large  n.  We  will  also  neglect  the  spontaneous  radiative  Auger  process  described  by  Aberg  *®, 
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which  is  not  expected  to  be  important  for  the  large  n values  of  interest  in  the  present  investi 


r 


gation.  The  photoexcitation  and  stimulated  emission  processes  which  can  occur  between  these 
states  will  be  systematically  treated  in  a subsequent  investigation.  Although  the  rate  for  radia- 
tionless capture  into  the  doubly-excited  state  a — /',  nl  decreases  as  n the  dominant  contri- 
butions to  the  dielectronic  recombination  rate  are  often  associated  with  large  values  of  n.  Au- 
toionization into  a low-lyini.  excited  state  i can  be  more  probable  than  autoionization  into  the 
ground  state  ig.  The  inclusion  of  this  additional  autoionization  process  has  been  found"  to 
result  in  a substantial  reduction  in  the  dielectronic  recombination  rates  for  certain  ionization 
states. 

The  stabilizing  transition  from  the  autoionizing  state  a to  the  state  J,  below  the  ionization 
threshold,  will  be  assumed  to  occur  only  by  means  of  the  spontaneous  radiative  decay  process 

;k'+(r)  (a)  _ + (j) 

with  a rate  denoted  by 

A^(J.a)  =A^(a^j). 

In  order  to  obtain  the  contributions  from  autoionizing  states  to  effective  ionization  rates,  it  will 
be  necessary  to  allow  for  the  inner-shell  collisional  excitation  process 

X + 0)  + e ■ — + (a)  + e- 

described  by  the  rate  coefficient 

Cg(a.j)  =C^(7-*a). 

It  has  been  argued  by  Bates  and  Dalgarno'^  and  by  Burgess  and  Summers^  that  the  collisional 
stabilization  rate  associated  with  the  inverse  of  the  excitation  process  (24)  can  become  compar- 
able with  the  spontaneous  radiative  decay  rate  only  at  very  high  densities,  because  the  inner- 
shell  de-excitation  process  involves  relatively  low-lying  levels.  It  is  also  apparent  for  the  same 


(22) 


(23) 


(24) 


(25) 


reason  that  the  corresponding  stimulated  radiative  decay  and  inner-shell  photoexcitation  pro- 
cess will  become  important  only  in  the  presence  of  intense  radiation  fields  at  relatively  short 
wavelengths.  These  additional  processes  will  be  treated  in  a future  investigation. 


The  total  collisional-radiative  transition  rate  IV (a.  a)  between  autoionizing  states  of  the 
z-times  ionized  ion  is  defined  by  an  expression  which  is  analogous  to  equation  (11). 

The  autoionizing  resonances  and  the  adjacent  continue  can  be  treated  as  separate  states 
provided  that  the  resonance  level  separations  are  large  compared  with  their  total  widths,  which 
are  determined  by  all  spontaneous  autoionization  and  radiative  decay  processes.  An  analysis  by 
Seaton  indicates  that  this  treatment  becomes  a better  approximation  with  increasing  z provid- 
ed that  the  radiative  decay  width  can  be  neglected.  Shore  has  shown  that  overlapping  reso- 
nances with  different  angular  momentum  quantum  numbers  can  still  be  treated  as  separate 
states. 


III.  EFFECTIVE  TRANSITION  RATES 


In  this  section  we  combine  and  generalize  previously  reported  statistical  theories  in 
order  to  derive  a set  of  rate  (time-dependent)  equations  for  the  population  densities  of  the 
low-lying  levels  which  contains  effective  transition  rates  describing  indirect  transitions  through 
the  more  highly-excited  bound  and  autoionizing  levels.  The  number  of  low-lying  levels  which 
are  to  be  explicitly  included  in  these  rate  equations  will  depend  on  the  problem  of  interest. 
Bates,  Kingston,  and  McWhirter^  allowed  for  the  time  dependences  of  only  the  ground  state 
population  densities  and  introduced  definitions  of  the  effective  recombination  and  ionization 
rates  through  the  bound  excited  states.  Burgess  and  Summers^  extended  these  definitions  by 
including  dielectronic  recombination  into  the  excited  levels  in  addition  to  direct  radiative  and 


three-body  recombination,  but  transitions  between  autoionizing  levels  where  neglected.  In  a 


between  a group  of  low-lying  levels,  which  have  relatively  long  relaxation  times;  but  processes 
involving  autoionizing  states  were  neglected. 

The  rate  equations  obtained  in  previous  investigations  describe  transitions  between  levels 
in  only  two  adjacent  ionization  stages.  By  taking  into  account  transitions  involving  bound  and 
autoionizing  levels  in  three  consecutive  ionization  stages,  we  are  able  to  derive  a more  general 
set  of  rate  equations  for  the  low-lying  level  populations  which  contains  three  different  general- 
ized effective  transition  rates.  The  familiar  corona  model  equations  for  the  ground  state  popu- 
lation densities  will  be  recovered  by  making  the  assumption  that  all  excited  states  decay  by 
only  spontaneous  radiative  or  autoionization  processes.  At  densities  for  which  only  the 
highly-excited  levels  important  in  dielectronic  recombination  are  appreciably  depopulated  by 
collisional  processes,  the  ionization  rate  remains  essentially  unaltered  but  the  effective  rate  of 
recombination  must  be  described  by  a collisional-dielectronic  recombination  rate  coefficient. 


Using  the  notation  introduced  in  the  preceding  section  for  the  various  transition  rates 
and  rate  coefficients,  the  rate  equations  satisfied  by  the  bound  level  population  densities  N(j) 
of  the  z-times  ionized  ion  can  be  written  in  the  form 


= - 1 00.;')  NW) 

J 

+ (j,  i)  N(i)  + Y.  ^0.  k)  Nik)  N^ 

I k 

+ L 0,  a)  N(a)  ij.  b)  Nib). 

“ » (26) 


where 

00.;')  - - Wij.j')  ior  j pi  j' 
and 


(27) 
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(28) 


QU.j)  - I WU'J)  + ^Jl5(/,y) 

j'^J  1 ' 

^a(k,J)  + ^Cg(a,J)  + C„p(A.  j) 


Equation  (26)  describes  the  population  of  the  bound  state  j due  to  bound-bound  transitions, 
direct  recombinations,  direct  ionizations,  stabilizing  radiative  transitions,  and  autoionization 
processes.  Equation  (28)  gives  the  total  depopulation  rate  Q{j,  J)  which  results  from  bound- 
bound  transitions,  direct  ionizations,  direct  recombinations,  inner-shell  collisional  excitations, 
and  radiationless  captures.  The  transitions,  involving  bound  and  autoionizing  states  in  the 
three  consecutive  ionization  stages,  which  are  taken  into  account  in  the  determination  of  N{j) 
are  illustrated  by  the  vertical  lines  in  Figure  1.  The  rates  or  rate  coefficients  associated  with 
the  various  excitation  and  de-excitation  processes  are  also  indicated.  Neutral  atoms  and 
completely-stripped  ions  clearly  represent  special  cases  for  which  only  two  adjacent  ionization 
stages  are  involved. 


The  autoionizing  level  population  densities  N{a)  of  the  z-times  ionized  ion  satisfy  the 
rate  equations 

dN(a) 


dt 


- ^ Q{a.  a')  N(a') 


£ i)  NO)  -I-  £ Cg(aj)  N(J) 


Ng. 


(29) 


where 


and 


Q(a,  a')  - — fV(a,  a')  for  a ^ a' 


(30) 


Qia,  a)  — £ IVia'.a)  + ^AgO.a)  -I-  ^A,(J.a). 


(31) 
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Equation  (29)  describes  the  population  of  the  autoionizing  level  a due  to  transitions  from  near- 
by autoionizing  levels,  radiationless  captures,  and  inner-shell  collisional  excitations.  Q (a,  a)  is 
the  total  depopulation  rate  arising  from  transitions  to  nearby  autoionizing  levels,  autoionization 
processes,  and  stabilizing  radiative  transitions. 


Following  the  procedure  of  Bates,  Kingston,  and  McWhirter^  we  shall  retain  the  time 
derivatives  for  only  a group  of  low-lying  bound  levels,  consisting  of  the  ground  level  and  a cer- 
tain number  of  metastable  levels,  which  we  shall  denote  by,  for  example,  j and  j'.  The  more 
highly-excited  bound  levels,  which  will  be  denoted  by  j and  j',  will  be  assumed  to  have  their 
steady-state  (statistical  equilibrium)  population  densities.  It  is  usually  assumed  that  ionization 
and  recombination  processes  occur  predominantly  from  the  low-lying  bound  levels,  because 
the  highly-excited  level  populations  are  negligible  in  comparison.  However,  in  order  to  obtain 
the  collisional-dielectronic  recombination  coefficients  associated  with  recombination  from  the 
r-times  ionized  ion,  it  will  be  necessary  to  allow  for  ionization  from  the  highly-excited  bound 
states  k produced  by  dielectronic  recombination. 


The  statistical  equilibrium  population  densities  N(J)eq  of  the  highly-excited  bound  levels 
J are  given,  in  terms  of  the  population  densities  of  the  low-lying  bound  levels  i,J,  and  k and 
the  autoionizing  levels  a and  b,  by 

-1  iT 


N(J)eq  = - I I C?  ■'  0./)  QiJ'J)  N{j) 

j'  i 


(jj')  « W,  i)  N(i)Ng 


+ IIQ-'  UJ')  Sir,  k)  N(k)  N, 


/'  * 


+ O'.y')  AfW.  o)  N(a) 


(32) 


In  the  next  section  we  will  describe  the  techniques  which  must  be  used  to  invert  the  rather 
large  matrix  Q,  which  is  defined  within  the  subspace  of  the  highly -excited  bound  levels  J. 
Analogous  expressions  are  obtained  for  N{i)g^  and  N(k)g^. 


We  will  assume  that  all  autoionizing  levels  a have  their  statistical  equilibrium  population 
densities  N(a)g^  given  by 

a f 

+ L L 0 (a.  a')  C,  (a'.J)  N(J)  A^. 

« j (33) 

When  transitions  between  the  autoionizing  levels  are  neglected,  the  first  term  on  the  right 

hand  side  reduces  to  the  usual  result  for  the  population  densities  of  the  doubly-excited  levels 

formed  by  radiationless  capture.  The  second  term,  which  gives  the  contributions  from  inner- 

shell  collisional  excitations,  is  important  only  for  relatively  small  values  of  the  outer  electron 

principal  quantum  number  n.  An  analogous  expression  is  obtained  for  Nib)^^  in  terms  of  the 

population  densities  N(j)  and  N(k). 


After  expressing  the  population  densities  of  the  more  highly-excited  bound  and  autoion- 
izing levels  in  terms  of  the  population  densities  of  the  low-lying  levels  in  three  consecutive 
ionization  stages,  we  obtain  a reduced  set  of  rate  equation  which  can  be  written  in  the  form 


dNU) 

dt 


- I N{J') 


+ ')  N(i)  Ng 


+ £ k)  Nik)  N,. 

k 


(34) 
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The  three  effective  transition  rates  which  appear  in  these  equations  are  rather  small  matrices 


which  are  defined  only  within  the  subspace  of  the  low-lying  levels  explicitly  included. 

The  matrix  elements  with  j ^ 7' are  the  effective  bound-bound  transition  rates 

between  the  low-lying  levels  of  the  ion  while  the  diagonal  elements  Qgffij.j)  represent 

the  total  effective  loss  rates  from  these  levels.  The  general  expression  obtained  for  QeffU.j') 
can  be  written  in  the  form 

QeffU.j)  - Ot-O./)  - II  OO.7O  Q ■'  (IT)  Qt(J'J') 

J T 

- I I fV,  SO,  k)  Q -1  (k.  k')  Wg  aj-ik'J') 

* *' 

- 1 1 « 0,  T)  Q (I  ?)  Ng  St  (Tf). 

i T (35) 

where 

Qt(J.J')  -Q(JJ')  -'LTAgUa)  Q-Ha.a')  NgCgia'J') 

» a 

- I I /fa  (J.  b)  Q "*  (b.  b')  Ng  C,,g(b'.r). 

* * (36) 

and  the  corresponding  definitions  of  ajik.j)  and  SjU,])  are  given  by  expressions  which  are 
analogous  to  equations  (38)  and  (40). 

The  total  transition  rate  QjU.j')  is  obtained  as  the  sum  of  the  direct  term  Q(j.j)  and 
two  additional  contributions  describing  indirect  transitions  through  the  autoionizing  levels  a 
and  b in  two  adjacent  ionization  stages.  The  third  term  in  equation  (36),  which  describes  au- 
toionization following  radiationless  capture,  reduces  to  the  result  obtained  by  Seaton  when 
transitions  between  the  autoionizing  levels  6 are  neglected.  The  second  term,  which  represents 
stabilizing  radiative  transitions  of  autoionizing  levels  a formed  by  inner-shell  collisional  excita- 
tion, has  not  been  discussed  previously. 
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J 


The  general  expression  for  which  includes  the  contributions  from  indirect 

transitions  through  the  highly-excited  bound  levels  in  three  consecutive  ionization  stages,  be- 
comes equivalent  to  the  definition  given  in  the  recent  paper  by  Burgess  and  Summers^  when 
only  the  direct  transition  rate  Q(J.J')  is  included  in  Qj-U.J')  and  the  last  two  of  the  three 
double  summations  are  omitted.  The  last  double  summation,  which  describes  recombination 
following  ionization,  is  not  expected  to  be  important  for  highly-excited  states  /.  However,  the 
preceding  double  summation,  which  describes  ionization  from  highly-excited  states  it  populated 
by  recombination,  will  be  important. 


The  general  expression  obtained  for  the  effective  recombination  coefficient 
be  written  in  the  form 

')  “ arU.  i)  “ £ Z Q UJ')  ocT^j' 

J T 

- I QtU'.O. 

' <■ 


(37) 


where 


“r  0-  ')  = “ (J.  ')  + Z Z O'  o)  Q ' (a.  a')  C.apia  . '). 

“ o (38) 

When  only  the  direct  recombination  coefficient  a(J,  i)  is  included  in  aj-ij,  /),  the  first  two 
terms  in  equation  (37)  correspond  to  the  generalized  collisional-radiative  recombination 
coefficient  recently  defined  by  Burgess  and  Summers.*  The  last  term,  which  describes  recom- 
bination from  highly-excited  levels  i produced  by  collisional  excitation,  will  not  usually  be  im- 
portant. The  definition  of  the  collisional-dielectronic  recombination  coefficient  between 
ground  states  which  was  introduced  in  the  earlier  paper  by  Burgess  and  Summers^  can  be 
recovered  by  neglecting  transitions  between  the  autoionizing  levels  a. 
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The  general  expression  obtained  for  the  elective  ionization  coefficient  is 


S.ffU  k)  - Sj-U.  k)  -El  00.7)  Q 0,7')  Sj-W.  k) 

' r 

“EE  50.  A)  Q a,  P)  QfiP.  k). 

* * (39) 


where 

Sj-U.  k)  = S(j.  k)  + E Z 0.  *)  Q ■'  ()>,  *')  (b'.  A). 

* ft  (40) 

This  effective  ionization  rate  coefficient  did  not  appear  in  previously  reported  statistical 
theories,  beca.  e of  the  inherent  restriction  to  only  two  consecutive  ionization  stages.  Gold- 
berg, Dupree,  and  Allen  have  found  that  the  contributions  from  autoionization  following 
inner-shell  collisional  excitation,  which  are  represented  by  the  second  term  in  equation  (40), 
can  be  comparable  to  the  direct  ionization  coefficient  S(j,  A). 

At  very  low  densities  where  practically  all  excited  states  decay  by  spontaneous  radiative 
or  autoionization  processes  in  times  that  are  short  compared  with  their  collision  times,  all  ioni- 
zation and  recombination  processes  can  be  assumed  to  originate  from  the  ground  states.  If  the 
ground  states  of  the  three  consecutive  ionization  stages  are  denoted  by  /^,  and  A,„  the 
effective  recombination  and  ionization  coefficients  apjjijg,  ij  and  A„)  which  are  ob- 

tained from  the  general  expressions  are  found  to  reduce  to  the  anticipated  corona  model 
definitions 


a,,  (z  -I-  1 — c)  = E 


« (7,  'o  ) + E 


AXi  «)  ',.) 


7 {a)  + A,(a) 


(41) 


and 


(z  — 1 — ' z) 


? Ajb)  + A,{b) 


(42) 


in  which  the  sum  over  j now  includes  ail  bound  states  of  the  z-times  ionized  ion. 
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If  use  is  made  of  the  identities 


1 -I 


(j.  b) 


A,  {k.  b) 


y A„{b)  A,(b)  ^ A„(b)  A,(b) 


(43) 


and 


1 -I 


Af  (J,  a) 


AaO’^) 


f A Ja)  + A,  (a)  ^ A„  (a)  + A,  (a)  ’ 

the  effective  loss  rate  QfffUoJo^  obtained  from  the  general  expression  can  be  reduced  to  the 


simple  result 


QeffOo-Jo^  “ [“c  (2  — Z - 1)  + Sj.  (z  — 2 + 1 ) ]. 


(45) 


which  is  the  sum  of  the  corona  model  recombination  and  ionization  rates  defined  by  analogy 
with  equations  (41)  and  (42). 

We  have  demonstrated  that,  in  the  corona  model  approximation,  the  recombination  and 
ionization  rates  which  appear  in  the  total  loss  rate  are  defined  in  the  same  manner  as  the  pro- 
duction rates  arising  from  transitions  originating  from  the  two  neighboring  ionization  stages.  If 
the  ground  state  population  densities  are  now  denoted  simply  by  N(z),  the  time-dependent 
equations  simplify  to  the  familiar  corona  model  relationships 
dN(z) 


dt 


— la^fz— *z  — 1)  -t-Sflz— z-l-l)l  N{z)  Ng 


+ a^{z  + \ ^ z)  N(z  + \ ) N,  -I-  (z  - 1 - z)  ^(z  - 1 ) N,. 


(46) 


The  requirement  that  ^ dN(z)ldi  — 0,  which  expresses  the  conservation  of  ions  under  the  as- 

Z 

sumption  that  the  excited  levels  have  negligible  populations  compared  with  the  ground  levels, 
is  most  naturally  satisfied  if  equation  (45)  is  valid. 

If  the  assumption  is  now  made  that  only  recombinations  can  produce  highly-excited 
states  which  can  be  appreciably  depopulated  by  collisional  processes,  the  effective  ionization 
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rate  coefficient  kg)  is  still  given  by  equation  (42).  However,  the  effective  recombina- 
tion rate  coefficient  afgUo-  <o^  8'ven  by  the  collisional-dielectronic  recombination 

coefficient  defined  by 

(z  -1- 1 — ' 2)  ’“aj(jg.ig)  - ^ Y.  OO0.7)  0 (7.y')“r^-' ' '0 

j*joT*j„  (47) 

where  the  coefficients  o j {J,  ig ) for  recombination  into  the  levels  J are  defined  by  equation 


The  simplification  of  the  general  expression  for  the  effective  loss  rate  QeffUo-Jo^ 
much  more  tedious.  We  anticipate  that  the  contributions  representing  ionization  of  the 
highly-excited  states  ^populated  by  recombination  can  be  combined  with  the  corona  recombi- 
nation terms  to  yield  the  collisional-dielectronic  recombination  coefficient  a^^iz  — » z — 1)  “ 
«c</  which  is  defined  by  an  expression  analogous  to  equation  (47).  This  simplification 

can  be  accomplished  by  utilizing  the  identities 

1 0-  b)  Q (b.  b') 

b j 

= r b)  (?-'(*.  b') 


1 - Y.  NgS(jg,k)  Q-Hk.k  ) 


- Y Q^kg.k)  Q-Hk.k  ). 

k*k. 


The  effective  loss  rate  QeffUo’Jo^  found  to  be  the  sum  of  the  corona  model  ionization 

rate  S^(z  — z + 1)  and  the  collisional-dielectronic  recombination  rate  a,,^(z  — z - 1). 
The  ground  state  population  densities  now  satisfy  the  equations  obtained  from  (46)  after  re- 
placing both  corona  model  recombination  coefficients  with  the  corresponding  collisional- 
dielectronic  recombination  coefficients.  In  the  following  section,  our  calculations  will  be 
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described  for  the  collisional-dielectronic  recombination  coefficient  defined  by  equation  (47). 


IV.  CALCULATIONS 

In  this  section  we  describe  the  details  and  present  some  results  of  our  calculations  for  the 
excited  level  populations  and  the  collisional-dielectronic  recombination  rates  of  multiply- 
charged  Fe  ions.  The  present  calculations  represent  an  extension  of  our  previously  reported 
corona  model  calculations"  to  the  intermediate  density  region,  where  the  low-lying  level  po- 
pulations are  still  close  to  their  corona  equilibrium  values  but  the  highly-excited  levels  impor- 
tant in  dielectronic  recombination  are  depopulated  by  both  collisional  and  radiative  processes. 


We  first  consider  the  populations  of  the  doubly-excited  autoionizing  levels  specified  by  a 
— /',  nl.  Since  inner-shell  collisional  excitation  can  be  neglected  for  large  n values,  we  can  as- 
sume that  these  levels  are  populated  predominantly  by  ihe  radiationless  capture  process  (19). 
The  radiation  less  capture  rate  coefficient  C,a,(/  — a)  can  be  expressed  in  terms  of  the  rate 
A„  (a  — /)  for  the  inverse  autoionization  process  by  means  of  the  detailed  balance  relationship 


C„,(,^a)  (£h/I<b  T,) 

|£(/)  - £(a) 


3/2 


X exp 


Ag  (a  — /). 


(50) 


The  statistical  weights  associated  with  the  energy  levels  £(/)  and  £(a)  are  denoted  by  g(i) 
and  g(a),  respectively,  and  Ef^  = 13.6  eV.  The  other  symbols  have  their  conventional  mean- 
ing. 


For  large  values  of  n and  r the  autoionization  rates  Ag  (a  — /)  can  be  obtained  by  using 
the  quantum  defect  theory  relationship  derived  by  Seaton  ",  which  may  be  expressed  in  the 
form 
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A^ia-'  i) 


2(r  -H)^  g(i)  EW)  -£(/) 
^ gia)  E„ 


X £ <r  (/,  t,  /, 


E,  ' t,  -0' 


(51) 


where  <r(i,  , «,  /,)|  denotes  the  threshold  value  of  the  partial  wave  cross  section 

€,  “0 

for  the  electron  impact  excitation  / — / of  the  (z  + l)-times  ionized  ion.  The  electron  impact 
excitation  cross  sections  used  in  our  calculations  have  been  obtained  by  the  distorted  wave 
method  described  by  Davis,  Kepple,  and  Blaha 


The  decay  rate  for  the  stabilizing  radiative  transition  (22)  is  obtained  by  making  the  cus- 
tomary approximation 


/4,(a— y)  -,4, (/'—/„), 

in  which  the  influence  of  the  outer  «/-electron  is  neglected. 


(52) 


The  autoionization  processes  and  stabilizing  radiative  transitions  which  are  taken  into  ac- 
count in  the  present  calculations  are  associated  with  single-electron  electric-dipole  de- 
excitations /'  — / of  the  recombining  ion.  The  radiative  decay  rates  corresponding  to  a 
An,  ^ 0 transition  increase  as  with  increasing  z.  However,  the  dielectronic  recombination 
of  the  Ee  — Ee  ions  is  accomplished  primarily  through  — 3p  transitions,  for  which 
the  radiative  decay  rates  increase  only  linearly  with  z.  Consequently,  the  stabilizing  radiative 
decay  process  becomes  competitive  with  autoionization  only  for  relatively  large  values  of  the 
outer  electron  principal  quantum  number  (n  ~ 100).  Finally,  we  note  that  the  3</—  3/j  colli- 
sional  stabilization  rate  becomes  comparable  to  -4,  (3</—  3p)  only  at  = 10*^  cm  Equa- 
tion (23)  is,  therefore,  valid  for  the  density  range  considered  in  this  investigation. 

The  rate  coefficient  describing  collisionally-induced  angular  momentum  mixing  of  the 
outer  /lAelectron  state  is  given  in  the  Bethe-approximation  by 
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(«  , 327rV2  \alE„ 

1)  z„  j— ^ 


3 

/> 

n 

4 

21  + 1 

z + 1 

(n^  - />  ) 


g// 


1/2 


2w 


fl  (>-). 


(53) 


where 


y 


2 

{£h 

2(z  + !)£« 

nig 

kgT 

(A£,2  + r/  + ). 


(54) 


The  charges  on  the  projectile  and  target  particles  are  denoted  by  Zg  and  z,  respectively,  and  fj. 
is  the  reduced  mass.  The  argument  y of  the  exponential  integral  £]  is  made  finite,  as  suggest- 
ed by  Griem  by  taking  into  account  the  plasma  frequency  the  /-sublevel  energy  separa- 
tion A£/,  and  the  width  T,  due  to  the  finite  natural  lifetime.  A similar  result  obtained  by 
Pengelly  and  Seaton  '*  does  not  include  the  /-sublevel  separation,  which  can  be  important  in  a 
multiply-charged  ion.  In  our  calculations,  the  /-sublevel  separations  were  estimated  by  taking 
into  account  the  long-range  monopole-quadrupole  interactions  and  the  fine  structure  splittings. 


The  normalized  population  densities  p{a)  ■=  iy(a)/N(ig)Ng  of  the  doubly-excited  levels 
a - / , nl,  which  are  formed  by  radiationless  capture  from  the  ground  state  /„,  are  given  by 


p(u)  = £(?-'(a.  a')  C„,(a  , /,,). 

a 


(55) 


With  account  taken  of  the  angular  momentum  redistribution  processes  (9)  and  (10),  the  non- 
zero elements  of  the  matrix  Q are  given  by 

-Q(a.a')  =IP(a,  a')  ^ N^Cginl'^  nl) 


+ /V^+  C^+  (.nl'  — nl) 


(56) 


and 
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(57) 


Q(a.  a)  — A^(a)  + A^(a)  + ^ a). 

a 


where  a differs  from  a only  in  the  outer  electron  angular  momentum  quantum  number  / = / 
± 1.  The  sum  over  a in  equation  (55)  is  actually  a sum  over  all  /'  = 0,  a — 1.  The  rale 
coefficient  describing  the  production  of  the  singly-excited  bound  states  j “ ig,  nl  by  dielectron ic 
recombination,  which  corresponds  to  the  second  term  in  quation  (38),  is  expressed  in  terms  of 
the  normalized  doubly-excited  level  population  densities  p (a)  by 


aj  0.  'o)  “ Z P 

a 


(58) 


Because  of  the  approximation  given  by  equation  (52)  for  A,,  (a  — J),  the  sum  over  a is  actually 
a sum  only  over  the  different  excited  states  /'. 


In  our  previous  investigation^,  the  effects  of  angular  momentum  redistribution  due  to 
protons  were  included  by  using  the  quasi-static  and  linear  Stark  effect  approximations.  The 
treatment  in  terms  of  proton  impacts  adopted  here,  which  is  essentially  the  same  as  that  of  Sea- 
ton and  Storey  is  expected  to  predict  less  angular  momentum  redistribution  for  a given  value 
of  n.  When  it  becomes  necessary  to  account  for  the  angular  momentum  redistribution  pro- 
duced by  other  impurity  ions,  the  quasi-static  may  be  the  more  appropriate  of  the  two  extreme 
approximations. 

With  increasing  n,  the  autoionization  rate  decreases  as  n whereas  the  rate  for  collision- 
ally  induced  angular  momentum  mxing  increases  as  n*.  However,  the  radiative  decay  rate 
given  by  equation  (52)  is  independent  of  n.  Consequently,  the  collisional  processes  are 
expected  to  be  capable  of  establishing  a Boltzmann  distribution  with  respect  to  / only  for  much 
higher  values  of  n then  in  the  case  of  the  singly-excited  bound  states,  for  which  the  radiative 
decay  rates  decrease  as  n Accordingly,  the  less  rapid  collisional  processes  connecting 
doubly-excited  levels  associated  with  different  /^values  can  probably  be  neglected. 
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The  effects  of  collisionally-induced  angular  momentum  redistribution  are  illustrated  for 
three  different  densities  in  Figures  2 and  3,  in  which  the  coefficients  obtained  from  equation 
(58)  for  the  dielectronic  recombination  of  Fe  are  shown  as  functions  of  / for  n — 20  and 
30,  respectively.  Both  the  electron  and  the  proton  temperatures  were  taken  to  equal  to  100  eV, 
which  is  close  to  the  maximum  abundance  temperature  predicted  by  our  previously  reported 
corona  equilibrium  calculations".  The  results  obtained  at  Ng  - 10'°cm  are  indistinguish- 
able from  the  values  predicted  without  the  inclusion  of  angular  momentum  redistribution.  For 
n - 30  (Figure  3),  the  results  obtained  at  Ng  - 10*®  cm  are  close  to  the  limiting  values 
predicted  on  the  basis  of  a Boltzmann  distribution  with  respect  to  /.  Results  similar  to  those 
presented  in  Figures  2 and  3 have  been  obtained  for  C by  Seaton  and  Storey 

The  coefficients  for  the  dielectronic  recombination  of  Fe  obtained  from  equation  (58) 
after  summation  with  respect  to  / are  shown  as  functions  of  n in  Figure  4.  As  anticipated,  an- 
gular momentum  redistribution  is  found  to  become  more  effective  with  increasing  Ng  and  n. 
The  amplification  in  the  dielectronic  recombination  coefficients  can  be  directly  attributed  to  an 
increase  in  the  populations  of  the  higher-/  doubly-excited  levels,  which  have  larger  statistical 
weights  and  practically  zero  probabilities  of  formation  directly  through  radiationless  capture 
and  decay  through  autoionization.  The  rate  of  photon  emission  through  stabilizing  radiative 
transitions,  which  is  obtained  after  summation  with  respect  to  n,  is  found  to  be  enhanced  by 
about  a factor  of  3 over  the  low-density  result  by  a density  of  10  ’®  cm  This  effect  may  be 
observable  in  the  emission  spectrum  near  the  resonance  line  of  the  recombining  ion. 

We  now  consider  the  populations  of  the  singly-excited  bound  levels  J - ig,  nl  produced 
by  recombinations.  According  to  equation  (38),  the  total  recombination  coefficient  aj-(j,  ig) 
into  level  j is  the  sum  of  the  total  direct  recombination  coefficient  a ij,  ig)  and  the  dielectronic 
recombination  coefficient  a^(J,  ig)  defined  by  equation  (58).  The  total  direct  recombination 
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coefficient  is  obtained  from  the  direct  radiative  recombination  coefficient  and  the  three-body 
recombination  coefficient  by  means  of  equation  (16).  The  direct  radiative  recombination 
coefficient  — 7 ) can  be  estimated  by  using  the  expression  presented  in  our  previous 

paper  which  takes  into  account  the  presence  of  equivalent  electrons  in  the  outermost 
subshell  of  the  recombining  ion.  The  three-body  recombination  coefficient  03  (/„  J)  is  ob- 
tained from  the  electron  impact  ionization  coefficient  S,  (7  — ) by  means  of  a detailed  bal- 

ance relationship  analogous  to  equation  (SO). 


The  electron  impact  ionization  rate  coefficient  obtained  in  the  Bethe  and  unit-Gaunt- 
factor  approximations  may  be  written  as 


C , 2«ir 

SeO  3^ 


h\  1 


if  I (2/  -1-  1 ) 


X |(£^/y  kgTf)^ 


(>)  -(-  y//i  4 - y 


exp ( -y) 


- iE„/kBT,)- 


li  X 1 
yln4 


-jy  ^exp(-y) 


~\y  ’ exp  ( -y)  -I-  £1  (y) 


(59) 


where 


y - [£(/„)  -E{j)\lkg  T,. 

(60) 

The  accuracy  of  this  expression  is  expected  to  improve  with  increasing  n. 

The  relative  importance  of  dielectronic  recombination,  direct  radiative  recombination,  and 

three-body  recombination  is  illustrated  in  Figure  5 by  showing,  'he  individual  contributions  to 
^ a j-  (7  /„)  as  functions  of  n for  the  recombination  of  Fe  at  - 100  eV  and  N^.  - 10 

I 
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Direct  radiative  recombination  can  be  seen  to  be  negligible,  while  three-body  recombination 
prevails  over  dielectronic  recombination  only  for  n > 40. 

In  order  to  obtain  the  population  densities  of  the  singly-excited  levels,  it  will  be  necessary 
to  estimate  the  bound-bound  transition  rates  which  enter  into  the  construction  of  the  matrix  Q. 
The  rate  coefficient  describing  ^ 0 electron  impact  excitations  is  given  in  the  Bethe  ap- 
proximation by 


Q (nl  - n'.  r 


/ ± 1) 


f32rr^/2 

1 >/3 

if 

where 


X f(nl  -*  n'  /')  y exp  ( -y) 


X (In  4 -I-  exp  (ylfj  (y)  1, 


(61) 


y - lEin'D  - EinDVkgT^ 

(62) 

and  finl-*  n'  I')  is  the  oscillator  strength  for  the  dipole  transition  which  is  approximated  in 
our  calculations  by  using  the  hydrogenic  value.  The  de-excitation  rate  is  obtained  from  the  ex- 
citation rate  by  means  of  a detailed  balance  relationship.  « — • n ± 1 transitions  are  found  to  be 
the  most  rapid  ^ 0 collisional  processes.  For  the  w-levels  considered,  the  n — • n + 1 excita- 
tion processes  are  found  to  be  more  probable  than  the  direct  ionization  process. 


In  previously  reported  calculations  for  the  collisional-dielectronic  recombination 
coefficients  it  has  been  assumed  that  the  collisional  processes  are  rapid  enough  to  establish  a 
Boltzmann  distribution  of  the  singly-exicted  bound  level  populations  with  respect  to  / for  all 
values  of  n.  This  assumption  becomes  questionable  in  multiply-charged  ions,  where  the  fine- 
structure  splitting  of  the  /-sublevels  scales  as  z*.  In  this  investigation,  we  take  into  account  the 
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rates  for  collisionally-induced  angular  momentum  mixing  of  the  singly-excited  bound  states  j by 
using  equation  (S4)  for  both  protons  and  electrons. 

The  radiative  decay  rates  A,(nl  — n',  I ± i)  are  approximated  by  using  the  hydrogenic 
values. 

The  normalized  population  densities  p(.j)  “ of  the  singly-excited  bound 

levels  / ~ /q,  n/,  produced  by  recombinations  from  the  ground  state  ig,  are  given  by 

P 0)  - 2^  0 0.7  ) “ r 0 . <0 )■ 

7— (63) 

The  matrix  elements  of  ©are  given  by 

- q(7J')  = iy(jj‘) 

= Ng  Cg  («'  ' nl)  + yV^+  (n'  /'— > nl)  6 (/?',  n) 

+ A^(n’  l'—>  nl)  «[£(»'/')  - E(nl)]  for  j ^ j' 

(64) 

and 

Q(JJ)  - Wg  ig)  + ^ ^ J • D ■ 

i*j  (65) 

The  collisional-dielectronic  recombination  coefficient  is  given,  in  terms  of  the  normalized 
singly-excited  bound  level  population  densities  p (7),  by 

+ X ^Oo.y)  p(./)- 

j*Jo  (66) 

With  increasing  n,  collisional  processes  will  tend  to  bring  the  singly-excited  level  popula- 
tion densities  closer  to  their  Boltzmann-Saha  equilibrium  values.  Accordingly,  it  is  customary 
to  express  the  normalized  population  densities  p (7  ) in  terms  of  the  departure  coefficients 
b(nl),  defined  by 
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p(7) 


X exp 


£(/„)  -£(;) 


(2/  + 1)  *(«/), 


(67) 


which  approach  unity  for  sufficiently  large  n.  The  matrix  Q must,  therefore,  be  constructed 
only  within  the  subspace  of  the  /(/-levels  for  which  b(nl)  differs  substantially  from  unity.  Since 
direct  inversion  of  this  matrix  is  still  impractical  for  the  large  //-values  of  importance,  we  have 
adopted  the  matrix  condensation  and  Lagrangian  interpolation  techniques  developed  by  Bur- 
gess and  Summers'**  in  order  to  reduce  the  matrix  to  a manageable  size.  Briefly,  a representa- 
tive set  of  III  levels  is  selected,  and  the  complete  set  of  level  population  densities  is  expressed 
in  terms  of  the  representative  values  by  means  of  the  appropriate  Lagrangian  interpolation 
coefficients.  Because  of  the  necessity  to  interpolate  with  respect  to  both  //  and  /,  the  present 
interpolation  scheme  is  greatly  increased  in  complexity  in  comparison  with  the  one  used  by 
Burgess  and  Summers^*,  who  assumed  that  the  /-sublevels  have  relative  populations  given  by 
the  Boltzmann  distribution  for  all  values  of  //. 


The  departure  coefficients  binl)  for  recombination  onto  Fe  are  shown  in  Figure  6 as 
functions  of  / for  //  = 10,  20,  and  30  at  /V,,  - 10 and  10**  cm  The  characteristic  overpo- 
pulation of  the  ///-levels  which  is  produced  by  dielectronic  recombination  can  be  clearly  seen. 
The  assumption  that  a Boltzmann  distribution  with  respect  to  / is  established  for  all  //-values, 
which  is  inherent  in  all  previously  reported  calculations  of  is  found  to  be  questionable. 
The  departure  coefficients  for  / — 0,  which  are  illustrated  as  functions  of  //  in  Figure  7,  show 
that  the  approach  to  Boltzmann-Saha  equilibrium  with  increasing  V,  and  //can  be  rather  slow. 

The  collisional-dielectronic  recombination  coefficient  for  the  recombination  of  Fe  is 
shown  in  Figure  8 as  functions  of  temperature  for  three  different  densities.  The  results  ob- 
tained for  — 10  '**  cm  are  essentially  the  same  as  our  previously  reported  " corona  model 
values,  while  the  results  obtained  for  A,,  - 10 cm  are  found  to  be  reduced  from  the  coro- 
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na  model  values  by  about  an  order  of  magnitude  due  to  the  effects  of  collisional  processes  on 
the  highly-excited  bound  levels  of  the  recombined  ion.  A similar  figure  has  been  presented 
for  Fe  by  Burgess  and  Summers^. 

The  collisional-dielectronic  recombination  rate  coefficients  obtained  for  the  recombination 
of  fe  — fe'*''^  at  S\.  — 10 cm  which  are  tabulated  in  Table  I as  functions  of  tem- 
perature, are  found  to  be  reduced  from  their  corona  model  values  by  at  least  an  order  of  mag- 
nitude. Recently*^,  the  effective  rate  coefficients  for  the  recombination  of  fe  — fe"*"” 
have  been  deduced  experimentally  by  analyzing  the  time-dependent  line  emission  spectra  pro- 
duced by  iron  ions  injected  into  a theta-pinch  plasma  with  an  electron  density  slightly  greater 
than  10*^  cm  Their  preliminary  analysis,  based  on  a time-dependent  corona  ionization- 
recombination  model,  indicates  that  the  effective  recombination  rate  coefficients  are  indeed 
lower  than  predicted  by  our  corona  model  calculations". 

V.  CONCLUSIONS 

In  this  investigation  we  have  generalized  previously  reported  statistical  theories  to  obtain 
the  contributions  from  highly-excited  bound  and  autoionizing  levels  to  effective  transition  rates 
between  the  low-lying  bound  levels  in  three  consecutive  ionization  stages.  We  have  also 
demonstrated  that  the  familiar  corona  model  equations  can  be  recovered  at  very  low-densities. 
Ai  densities  for  which  only  the  highly-excited  levels  involved  in  dielectronic  recombination  are 
appreciably  depopulated  by  collisional  processes,  the  effective  ionization  rate  is  practically  unal- 
tered but  the  effective  recombination  rate  can  be  substantially  reduced.  Calculations  have  been 
carried  out  for  the  excited  level  populations  and  the  collisional-dielectronic  recombination 
coefficients  of  Fe  - Fe  ions,  which  recombine  through  very  large  values  of  the  outer 
electron  principal  quantum  number.  To  the  best  of  our  knowledge,  the  present  calculations  are 
the  first  in  which  explicit  account  has  been  taken  of  the  transition  rates  for  collisionally- 
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induced  angular  momentum  mixing  of  both  the  singly-excited  bound  states  and  the  doubly- 
excited  autoionizing  states.  In  a future  investigation,  more  extensive  results  for  collisional- 
dielectronics  recombination  coefficients  will  be  presented  together  with  new  results  for  the  ioni- 
zation equilibrium  and  the  radiative  energy  loss  rates  for  various  low-  and  high-Z  elements. 
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Table  1 — a (cm®  sec"^),  = 10^®  cm“®. 


^oi\oTt{eV) 

F«*8 

Fe*^ 

fe+lO 

life+13 

1.6 

0.61(-10) 

0.48(-10) 

0.47(-10) 

0.10(-9) 

0.11(-9) 

0.13(-9) 

1.8 

0.28(-10) 

0.21(-10) 

0.27(-10) 

0.51(-10) 

0.40(-10) 

0.27(-10) 

2.0 

O.ll(-lO) 

0.16(-10) 

0.17(-10) 

0.24(-10) 

0.18(-10) 

0.13(-10) 

2.2 

0.12(-10) 

0.97(-ll) 

0.89(-ll) 

0.91(-11) 

0.75(-ll) 

0.67(-ll) 

2.4 

0.2U-10) 

0.63(-ll) 

0.54(-ll) 

0.53(-ll) 

0.41(-11) 

0.29(-ll) 

2.6 

0.22(-ll) 

0.17(-11) 

0.18(-11) 

0.34(-ll) 

0.24(-ll) 

0.17(-11) 

2.8 

0.17(-11) 

0.12(-H) 

0.12(-11) 

0.16(-11) 

0.13(-11) 

O.ll(-ll) 

3.0 

0.14(-11) 

0.72(-12) 

0.71(-12) 

0.92(-12) 

0.78(-12) 

0.62(-12) 

31 


Fig.  2 - The  rate  coefTicientt  for  the  dielectronic  recombination 
of  into  the  nWeveU  of  Ft**  corresponding  to  ft  = 20. 
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